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ABSTRACT 

0^ 

An analytic expression is given for the growth function for linear perturbations in a 
Qh' low-density universe made flat by a cosmological constant. The result involves elliptic 

. integrals but is otherwise straight-forward. 

(N 

' Subject headings: cosmology: theory - large-scale structure of the universe 

> 

ly-^ I In linear cosmological perturbation theory, the evolution of perturbations in the absence of 

' pressure is reduced to the superposition of two modes with fixed time dependence and arbitrary 

■ spatial dependence ( Peebles| 1980). The behavior of these modes in time is a function of the 



. background cosmological parameters. The functions of time that describe this behavior are called 

. the growth functions and denoted Di and D2- 

^ ■ 

As shown by pieath] (1977), the growth functions in cosmologies where pressure fluctuations 



O I are negligible can be written as 

S2 ■ - ^ . , . r da' 



Di{a) oc H{a) / , (1) 
D2{a) oc H{a), (2) 



CI ■ H{a) = V^oa-'^ + ^Ra-^ + ^A- (3) 

Here, a is the expansion scale factor of the universe, chosen so that a = 1 today; this will be our 
time coordinate. The redshift z is then a~^ — 1. is the density of non-relativistic matter in the 
universe today in units of the critical density. J7a = ^/"^Hq is the cosmological constant A relative 
to the present-day Hubble constant Hq. Qji is the curvature term and is equal to 1 — 0,q — JIa- We 
assume that relativistic matter is a negligible contributor to the density of the universe. Finally, 
H{a) is proportional to the Hubble constant at epoch a, scaled so that H = 1 today. We are 
interested in the behavior of Di{a), as this is the mode whose amplitude grows in time. However, 
the overall normalization of Di is merely a convention. 



The analytic form for Di in the A = case ( IWeinbergI 1972; |Groth fc Peelolei 1975; [Edwards fc| 



Heath| 1976) is widely known, but no analytic solution for the case of A 7^ has been presented in 



the literature. Instead, workers integrate equation (|^) numerically or use the approximations given 
by [Lahav et al.| (1991) or Carroll et al.| (1992). However, the integral can be done in terms of elliptic 
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integrals, which of course are implemented as special functions in many numerical packages. Here, 
we restrict ourselves to the case of a flat, low-density universe, i.e. = 0, f^o < 1, = 1 — ^o- 
This is the case commonly used in practice and the solution can be stated rather cleanly. 



We adopt a normalization for Di as 



Di{a) 



H(a) 



da' 



a'^H{a') 



(4) 



This is chosen so that Di{a) a as a —> Q. Defining v = a ^Oo/(1 — ^o) and manipulating the 
integral yields 



Di{a) = -av^Jl + v^ / — 



u du 



+ n3)3/2 

d 

_.Vl + t;3 lim lim(-2)— / 

2 voa^cob^i do Jv 

5av\/ 1 + lim 



-Of A 



u du 
Vb + u 



u du 



+ 



(5) 
(6) 
(7) 



Doing this last integral (e.g. using 3.139.8 in [Gradshteyn fc Ryzhik| 1994) yields the final answer 

Di{a) = a X di i 



where 



diiv) 



3 

Arccos 



V + l-y/3 



E{P,sm 75°) 



1 



3 + ^3 



F(/?,sin 75°; 



+ 



+ l + \/3' 



(8) 

(9) 
(10) 



Here, F and E are incomplete elliptic integrals of the first and second kind, using the definitions 
F((/., k) = /(f (1 - sin^ Oy^l'^de and E{(j3, k) = (1 - A;^ sin^ Of/'^de; note that conventions vary 
on whether k is squared or not. Also note that (3 can exceed vr/2; some numerical implementations 
of the elliptic integrals do not handle this properly. In this case one can use the identity F{(j), k) = 
2F{tt/2, k) — F{tt — (f), k) and likewise for E to compute for (3 > Tr/2. 

As a — > 0, u ^ +00 and di ^ 1 — (2/11)^^^ + (16/187)f ^® + 0{v~^). This is not obvious from 
equation (|9|); the terms of order v"^ and v cancel. This suggests that for numerical evaluation for 

V ^ 1, one should be watchful for loss of precision. Alternatively, for f ^ 5 one can use the large 

V expansion given above. As o — > oo, v ^ and Di = 1.4373 -^170/(1 — ^^o)- 



The approximation of Carroll et al.| (1992, adapted from Lahav et al.| 1991) that 



nl''^ - r?A + (1 + l^o/2)(l + 1^a/70) 



(11) 



- 3 - 



can be compared to di{v) using Q\ = 1 — = {v^ + The fit has fractional errors better than 
2% over the range z > and > 0.1 (i.e. v > 0.48) but deviates for lower Qq. A plot of di{v) 
and the Carroll et al.| (1992) fit is shown in Figure 1. The fit has similar accuracy for the A = 
case even to much lower r^Q. 



For completeness, the solution for the expanding epoch of a flat universe with l^o > 1 is 
Di{a) 



a X di i - ' 



fin 



diiv) 



7 



1 



V3 



F(7,sinl5°) -£(7, sin 15°) 



+ 



Arccos ■ 



1 - Vs 



u = 1 is the time at which the universe reaches maximum expansion. As v 00, di 
(2/ll)v-3 + (16/187)^^6 ^ o{v-^). 



(12) 
(13) 
(14) 
1 + 



As a convenience to the reader, we collect other exact equations relevant to low-density, flat 
universes from the literature ( [Edwardl 1972, 1973; Dabrowski &: Stelmach 1986, 1987; [Weinberg 
1972, 1989). The time-redshift relation may be written as 



t{z) 



1 



sinh"^ {{1 + z) 



-3/2, 



'1 - Qr 



fin 



(15) 



The conformal time r] = J a ^dt is 

r]{z) = 



J_ fC^+^r' da 
Ho Jo a^ia) 
1 



o 1 -1/6 

nlil-flo) F[/3(z),sin75°] 



(3{z) 



Arccos 



l + z + (l- V3)^/j)o^ -1 



l + z + (l + V3)V^o^ -1 



(16) 
(17) 

(18) 



where H{a) is defined in equation (|3|). The coordinate radius out to redshift z is ri{z) = c[r]{0) — 
r]{z)], where c is the speed of light and the flatness of the universe has been used. Then the 
luminosity distance is (1 + z)ri{z), the angular diameter distance dA to an object of given 
physical size is (1 + z)^^ri{z), and the comoving volume per unit redshift per unit solid angle is 



dVc/ dz dJl 



rliz) 



Ho ^l-fto + noil + z)-^' 



(19) 



A quantity related to the growth function is / = {a/ Di){dDi/da), used to relate density 
perturbations to velocity perturbations. Assuming the normalization in equation (^, one finds 
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(peeble^ 1984) 

s _ ^0 / 5a _ 3 

~ {l-no)a^ + no\2Di{a) 2 

In summary, the expression presented here for the growth function is marginahy more com- 
phcated than the A = case, but if one evaluates the elhptic integrals using available numerical 
packages, the exact function is quite tractable. 
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Equation (9) 
Carroll et al. Fit 



.1 



1 

V 



Fig. 1. — The function di{v) from equation (^) {solid line). The growth function Di{a) is 
adi{a~^ y^Tlo/JT^^TlQ). The approximation of Par roll et al. (1992) (dashed line), using 1 — ^Iq = 
(v^ + fits quite well in the observationally relevant portion of parameter space. 



